Homogenization of linear transport equations in a stationary 

ergodic setting 

Anne-Laure Dalibard* 



Abstract 

We study the homogenization of a linear kinetic equation which models the evolution of the 
density of charged particles submitted to a highly oscillating electric field. The electric field and the 
initial density are assumed to be random and stationary. We identify the asymptotic microscopic 
and macroscopic profiles of the density, and we derive formulas for these profiles when the space 
dimension is equal to one. 

1 Introduction 

This note is concerned with the homogenization of a Unear transport equation in a stationary ergodic 
setting. The equation studied here describes the evolution of the density of charged particles in a rapidly 
oscillating random electric potential. This equation can be derived by passing to the semi-classical limit 
in the Schrodinger equation (see [9], [TT], and the presentation in [7]). Our work generalizes a result of 
E. Frenod and K. Hamdache (see |7j) which was obtained in a periodic setting. The strategy of proof we 
have chosen here is different from the one of |7j , and allows us to retrieve some of the results in |3 in a 
rather simple and explicit fashion. 

Let us mention a few related works on the homogenization of linear transport equations; we emphasize 
that this list is by no means exhaustive. In p], Y. Amirat, K. Hamdache and A. Ziani study the 
homogenization of a linear transport equation in a periodic setting and give an application to a model 
describing a multidimensional miscible flow in a porous media. In [3] (see also f8]), Laurent Dumas and 
Frangois Golse focus on the homogenization of linear transport equations with absorption and scattering 
terms, in periodic and stationary ergodic settings. And in [4], Weinan E derives strong convergence 
results for the homogenization of Hnear and nonlinear transport equations with oscillatory incompressible 
velocity flelds in a periodic setting. 

Let us now present the context we will be working in : let (fl,T,P) be a probability space, and let 
{'''x)xeR'^ be a group transformation acting on fl. We assume that preserves the probability measure 
P for all X G K^, and the group transformation is ergodic, which means 

VA e T, {t.^A ^AVxeR^^ P{A) = or l) . 

The periodic setting can be embedded the stationary ergodic setting (see [13]). We will denote by E[-] 
the expectation with respect to the probability measure P; in the periodic case, we will write (/) rather 
than E[f] to refer to the average of / over one period. 

We consider a potential function u — u{y,uj) E L°°(IR^ x fi) which is assumed to be stationary, i.e. 

u{y + z,uj) = u{y, t^oj) V(y, w) G x x 

Moreover, we assume that < u{y,uj) < Umax = supu for all y G M^,cj e fl, and ■ti(-,w) g W^^^^{M^) 
for almost every u G Q,, so that \'yu{y,Lu) is well-deflned and locally Lipschitz continuous with respect 
to its y variable. 
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Let — f^{t, X, £,,uj), (t > 0, X E M.^ , ^ G R^, uj E il) he the solution of the transport equation 

r dtnt,x,^,to)+^-VJ%t,x,t^)~^^yu{l,u;) ■V^r{t,x,i,u)^0, 
\ r(t = 0,x,C,c.) = /o(a:,f,e,u;). 

Here, we assume that the initial data /o = /o(a;,y,C,t^) belongs to Ll^,{^^ x ]R|',L°°(K^ x ft)) and 
is stationary in y, i.e. 

/o(a:,y + z,<^,cj) = /o(a;, y, t^w) for all (x, y, z, ^, cj) e R"*^ x n. 

It is well-known from the classical theory of linear transport equations that for every uj E fl, there 
exists a unique solution of ^ in L^^{{0, oo), ^[^^(.(M;^ x M^)). The goal of this paper is to study the 
asymptotic behavior of as £ ^ 0. Thus, following [7], we define the constraint space K : 

Definition 1.1. Let 

e • V^/(y, u) - Vyu{y, to) ■ V^f{y, ^,uj) = (2) 
be the constraint equation, and let 

K:={f E Ll^iRf X R^,L\n)); f satisfies in V'{M.^ x ) a.s. in cu}. 

We also define the projection P onto the constraint space K, characterised by P{f) E K for f E 
^oci^i X R^,L^{n)) stationary, and 

I iP{f)^f)iy,^,Lo)g{y,C,Lo)d^dP{L,) = fora.e.yER^ 

for all stationary functions g E L°°{Ry x R^ x il) n K, with compact support in ^. 
(A more precise definition of the projection P will be given in the second section). 
Finally, we define as 

:= {/ E L]JM.^ X R^,L\n)y,3g E x R^,L\n)), f = P(g) -5}. 

Remark 1.1. Let us indicate that the constraint equation can easily be derived thanks to a formal two- 
scale Ansatz : indeed, assume that 

f^{t,x,^,u) K f (t,x,^,t^) ase^O; 

inserting this asymptotic expansion in equation we see that f necessarily satisfies the constraint 
equation 

Remark 1.2. Let f,g E L°° (Ry , L'^ (R^ xQ)) be stationary, and assume that f eK and g eK-^. Then 
for a.e. y E R^ , 

f f{y,^,u)giy,^,u)d^dP{uj)^0. 

This is a characterization o/K^ for the class of stationary functions in L°° [Ry , L'^ (R^ x Qj). 

Here, we provide another proof for the result of E. Frenod and K. Hamdache in |[7] in the "non- 
perturbed case". Our proof is based on the use of the ergodic theorem, and gives a more concrete insight 
of the projection P and of the microscopic behavior of the sequence Moreover, it allows us to retrieve 
the explicit formulas of the integrable case. 

The first result we prove in this paper is the following 

Theorem 1. Let fa E L^^iR^ x R^ x R^;L^{n)) stationary. 

Let f^ = f^{t,x,^,uj) be the solution of ([l]). Then for all e > 0, there exist f = f(t,x,y,S,,u!) and 
g = g{t,x;T,y,S^,u!), both stationary in y, and a sequence {r^(i, x, ^, a;)}e>o such that 

f''{t,x,£,,u}) ^ f (^t,x,^,S,,iS^ +g (^,x;^,^,^,u?j + r''{t,x,£,,uj) 

and : 



2 



• lk''llLj^^((0,oo)xR«xRf ,Li(f2)) ^0 fts e -> 0; 

• / e L^^((0,oo);ii^^(Rf xRf X R^;L^{n))), and f{t,x) £ IK for a.e. t>0,xe M^; 

• For all T>0, for all compact K cR^ xR^ xR^, 

sup WgWL^Kxn) < oo. 

0<t<T,0<T<T 

Moreover, g{t,x;T,-) e K-'- for a.e. {t,x,T) e (0, oo) x x (0,oo); 

• Microscopic evolution equation for g : for a.e. t,x (0,oo) x R^ , g{t,x; •) is a solution of 

1^ + e ■ Vy.g - Vyu ■ V^g = 0. (3) 

Moreover, for all T > 

as £ ^ 0. 

Li^^(R^xR|',Li(0)) 

• Macroscopic evolution equation : f and g satisfy 

9t(^)+e«(2/,C,^)-V.(^^)=0, (4) 

where 

:=P(0(2/,C,c.); 

• Initial data : 

f(t = Q,x,y,^,uj) = P{fo){x,y,£„uj), 
g{t = 0, x; r = 0, y, w) = [/o - P{/o)] {x, y, w). 

Before going any further, we wish to make a few comments on the above results. First, let us stress 
that it is not obvious that the function g is well-defined : indeed, let S{t) (t > 0) denote the semi-group 
associated to the macroscopic evolution equation and let T(t) (r > 0) be the semi-group associated 
to the microscopic evolution equation ^ . Then g is well defined if and only if, for all stationary function 
go = go{x,y,^,uj), for all t,T>0, 

T{t) [S{t)g^] = Sit) [T{T)go] . 

This identity follows from the fact that the speed S^^y,^,^) appearing in equation |[4]) is a stationary 
solution of (lU by definition of the projection P, and is thus invariant by the semi-group T(t). 

Next, let us explain briefiy the meaning of theorem [TJ The idea is the following : write /o as 
ft) = /o|| + /o-L, with /o||(a;, •) S IK and fo±{x, ■) G IK^ a.e. Then f^ can be written as /| + /j^, where 
/| (resp. ff) is the solution of equation ^ with initial data /o|| (x, |,^,a;) (resp. /o_l {x, j,£,,i^)). 
Theorem [T] states that 

fl-f(t,X,^,^,Lu)^0 

strongly in Lj^^ norm. In particular, there are no microscopic oscillations in time in this part of f^. We 
wish to emphasize that this result appears to us to be new. 

We now focus on the other part, namely An easy consequence of the theorem is 

/ fi{t,x,tuj)dt^O 
Jo 
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Ll^{R^;L^{n))) and for all T > 0. However, it would be wrong to think that /j^ vanishes 



in Ll^{{0,oo) X X R^,L^{n)), for instance. Indeed 
in Lj'gj,, and 



t X 

e e 



\\g{t = 0,x;T,?/)||ii(R«>,n) = ll/o±(a;, y)||Li(R«xn) 

as soon as fo±(x,y) S L^(R^ x fl) for almost every x,y. Consequently, if /o_l 7^ 0, then for all T > 
and for all compact K C K^, there exists a constant C > depending only on K, ||/oj_||li , and T such 
that ^ 

t X 



dt > c. 



Li(i<'xR?'xn) 



Hence does not vanish strongly in general. In other words, there are fast oscillations in time, due to 
the ill-preparedness of the initial data (i.e. fo{x, •) ^ K), but these oscillations do not cancel out as e 
vanishes. 

Let us now explain briefly here how our strategy of proof differs from the one of E. Frenod and K. 
Hamdache. The key of our analysis lies in the study of the behavior as e ^ of the Hamiltonian system 

Y'it,x,^,oj) = ~E'{t,x,^,iu), t>0 
EHt, y, e, to) = ^Vyu{Y^{t, x,tco),Lo), t>0 
Y^{t = 0, X, Lu) = X, E^{t = 0, X, oj) = {x, ^, iv) e M^^^ x n. 



fit, X, u;) = fo[ Y^t, X, C, c.), ^^iMil^, S-(t, y, u;),u;], 



Indeed, 



so that we can deduce the asymptotic behavior of from the one of (F^, ^^). And it is easily checked 
that 

y-(t,a;,^,a;) =er(|,f,e,c^) 



S-(t,y,^,a;)=S(|,f,^,a;j, 
where (Y, S) is the solution of the system 

Y{t,y,^,uj) = -E{t,y,^,uj), t>0 
^{t, y, uj) = Vj;u(y(t, y, w), w), t>0 

Y{t = 0, y, C, u;) = y, E{t = 0, zj, ^ u^) - t (y, ^, ^) G x ^. 



(5) 



Hence, in order to study the limit of as e — > 0, we have to investigate the long time behavior of the 
system {Y, S), and this will be achieved with the help of the ergodic theorem in the second section. 

In the case when iV = 1, we can give explicit formulas for ^}[y^ f , w); the proof of this formula in the 
stationary ergodic case is strongly linked to methods from the Aubry-Mather theory (see [5], [6], [12]), 
and thus also to the homogenization of Hamilton- Jacobi equations. Let us first recall the definition of 
the homogenized Hamiltonian H (see [lOj) 




A 



if IpI < E 
if IpI > E 



V2(Wmax - U) 



U) 



, where |p| = E -\/2(i 



u) + \ 



Proposition 1.1. Assume that N — I. 

Let (y, ^, w) g M X R X £7 such that H{y,^,uj) > u^ax- Assume that for all Q e M, for all w G il, 
there exists a Lipschitz continuous function ti(-, w), viscosity solution of 



H{y,Q + Vyviy,Lu),Lo)^H{Q) 
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such that 

^ _^ ^0 as |y| -> oo (6) 

U.S. in uj. 

Let P = P{y,^,Lu) G R such that H{P) = H{y,^,Lu) and sgn{P) = sgn{£,). Then 

e{y,^,u)^H'{P) 
Moreover, if L is the dual function of H , i.e. 

L{jJ,P,u)) = sup(p^ - H{y,^uj)) = hp\'^ - u{y,ui), 

and L is the homogenized Lagrangian, then 

P(L)(y,^,c.) = L(C»(2/,e,^).) 



In the periodic case, we will give another proof of the above result; the strategy chosen in that case 
is inspired from techniques and calculations in classical mechanics. It also allows to give a formula for 
for low energies in the periodic setting only: 

Proposition 1.2. Assume that N = I and that the environment is periodic. 
Let e R2 such that H{y,^) < u^ax- Then = 0. 

The organisation of this note is the following : in the second section, we derive some preliminary 
results on the long-time behavior of the system {Y, E) thanks to the ergodic theorem. Those will be 
useful in the proof of theorem [TJ to which is devoted the third section. Eventually, the fourth and last 
section is concerned with results in the integrable case, both in the periodic and the stationary ergodic 
settings. 

2 Preliminaries 

This section is largely devoted to the study of the long-time behavior of the Hamiltonian system (Y, S) 
defined by ([5]). First, notice that the Hamiltonian 77(y,^,w) :— ^ICP + u{y,Lj) is constant along the 
curves of the system (F,5), and if / € L°°{n,C^(R^ x Rf)) is stationary, then 

feK ^ f{Y{t,y,t^),my,^,^),^) = f{y,i:^) V(2/,C,c.) e x x l]. 
Indeed, for all / G L°°{Q,C\R^ x Mf )), we have 

^/ iY{t, y, e, Lo),Eit, y, ^) = {H, /} iY{t, y, C, c^), S(t, y, uj),lo), 
where {H, /} denotes the Poisson bracket of / and H, i.e. 

{H, /} {y, e, u;) - C • V,/(2/, ou) - Wyu{y, ^, uj) ■ W^fiy, ^, uj). 

Let us mention an easily checked property of the trajectories (Y, E) which will be used extensively in 
the rest of the article : for all (y, z, ^) G M^^, for all lu e fl, t > 0, 

Y{t, y, Tzuj) + z = Y(t, y + z,^,uj), 

E{t,y,^,T,uj) ^E{t,y + z,^,u;). (7) 

In the periodic case, this invariance entails that the hamiltonian system (F, S) can be considered as a 
dynamical system on the N dimensional torus [0, 27r)^. In this periodic setting, it is somewhat natural 
to introduce the semi-group of transformations {Tt)t>a on [0, 27r)^ x R^ given by 

Tt{y, = {Y{t, y, 0, y, 0), y G [0, 2^)^, ^ G R^. 
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According to Liouville's theorem, this semi-group preserves the Lebesgue measure on [0, 27r)^ x R^; 
moreover, we can construct a family of finite invariant measures on [0, 2tt)'^ x by setting mdy, ^) = 
'^H{y4}<c dy for c > (remember that the Hamiltonian is constant along the hamiltonian curves). 
This construction is the root of the ergodic theorem (see corollary 12. ip , and thus of the study of the 
long-time behavior of the system {Y,E). 

In the stationary ergodic setting, this construction can be generalized as follows : we define the 
transformation Tt : x Q ^ Rf x Q hy 

together with the family of measures 

:= 1-H(c,i^)<c dP{uj) 

where 7i(C, ui) := -t- u(0, uj). It is obvious that for all c € (0, oo), fic is a finite measure on R^ x fl. 

Notice that the "good" generalization to the stationary ergodic setting of the semi-group {%) is a semi- 
group which acts onR^ xfl rather than R^ xRf. Thanks to the group of transformations (Tx)xeR'^ ■> 
the transformations in can result in transformations in Ry , but the definition chosen here ahows us 
to define a family of finite invariant measures, whereas such a construction is rather difficult if one tries 
to define a semi-group acting on x R^ . This will be fundamental in the rest of the proof. 

Lemma 2.1. {Tt)t>o is a semi-group on R^ x il and preserves the family of measures fic- 

Proof. Let us first prove the semi-group property : let t,s E [0, oo), and {£,,uj) G R^ x il; then 

TtoTsi^u;) = Tt(S(s,0,^,c^),Ty(,^o,«,^)'^) 

= (S(i, 0, E{s, 0, £,, uj), ry(s^o,e,c^)^), ^0 

and using the properties ([7]) we deduce 

E{t,0,E{s,0,^,iu),TY(^sM.i,^)^) - S(t,r(s,0,e,c^),S(s,0,e,c^),w), 

= S(t + s,0,^,w) 



and 



'''Y{t,Y{s,0,^,uj),S(s,0,^,uj),uj)'^ 
'''Y{t+s,0,i,uj)'^ 



Thus 



Since it is obvious that Tq = Id, {Tt)t>o is a semi-group on x fi. 

We now have to check the invariance property; let F e L^{R^ x f2; fic) arbitrary. We set f{y, uj) 
F{^,TyU;) for (y,^,^) € x x fi, and we compute 



/ FiTt{^,uJ))d^l,{^,u;) = E 

JB"xO 



f{Y{t, 0, l^),^it, 0, ^, Uj),Uj)lH{Y{t,0,i,Lj),Sit,Q,i,uj),uj)<c d^ 



Since the probability measure P is invariant by the group of transformation Ty , and 

f{Y{t, y, ^), 2(i, y, uj), uj) = f{Y{t, 0, ^, TyUj), E{t, 0, ^, TyUj),TyUj), 
we have, for all y £ 

E [f{Y{t, 0, uj), E{t, 0, uj), uj)lH(Y{t,o,i,i.j),s{t,o,i,<.j),i.j)<c] = 

= E [f{Y{t, y, uj), E{t, y, ^, uj), w)lff(F(t,y,^,L^),H(t,y,5,<^),<^)<c] 
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Take an arbitrary function </) € L^{M.y), and write 



E 



dyd^cl,{y)f{Y{t,y,^,uj),E{t,y,^ 



We change variables in the integral in (y,^ by setting {x,v) — (Y{t,y,^,uj),'B{t,y,^,uj)); according to 
Liouville's theorem, the jacobian of this change of variables is equal to 1, and 

{x,v) ^ (Y{t,y,^,uj),E{t,y,^,uj)) ■^=^ {y,^) ^ {X{t,x,v,u;),V{t,x,v,uj)), 

where {X, V) is a solution of the Hamiltonian system 

V = ~Vu{X,Lu), 
{X, V)it^O,x,v)^{x,v). 

Observe that in the present case, we have simply 

X{t, X, V, Lu) = Y{t, X, ~V, Uj), 

so that {X, V) satisfies relations ^ . 
Hence 

dy d^(l){y)f{Y{t, y, ^, Uj),E{t, y, ^, Uj),Uj)lH(Y(t,y,i,u:),Eit,y,i,uj),uj)<c 
dx dv(p{X{t, X, V, Uj))f{x, V, Uj)1h(x,v,uj)<c 
dx dv(t){X{t, 0, V, T^Uj) + X)F{V, T.^Uj)lH{v,r^i^)<c 



so that 



R2N 



= E 
= E 
= E 
= E 

since the integral of (p is equal to 1 



[ F(Tt(e,c^))dMc(C,c^) 

dx dv(l){X{t, 0, V, TxUj) + x)F{v, TxUj)Ih{v,t,u)<c 

dx dv(l){X{t, 0, V, uj) + x)F{v, uj)lH{v,i^)< 
dv ( / (j){X{t,0,v,uj) +x) dx] F{v,uj)l-H{v,uj) 
dvF{v,uj)l-H{v,L,)<c = I F duc 



The following corollary is an immediate consequence of Birkhoff's ergodic theorem: 



□ 



Corollary 2.1. Let F e L^(M^ x il;^c)- There exists a function F € L^(R^ x fi; ^c) such that as 

T oo, 

F{Tt{tu^))dt^Fi^,L,) 
a.e. on x and in L (fic)- Moreover, F is invariant by Tf for all t > 0, and 



R" xO 



F dfic = 



F dfic- 



(8) 
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Additionally, if f = f{y, ^, uj) is the stationary function associated to F, that is, f{y, ^, lo) — TyUj), 
then f is invariant by the hamiltonian flow (F, S); precisely, for all {y,^,uj) S x H,, t > 

f{Y{t,y,^,Lo),E{t,y,^,co),Lo) = f{y,L^)- 



Proof. We only have to prove the invariance of / by the Hamiltonian flow; first, for y = 0, we have 

fiY{t, 0, t w), 0, uj),Lu) = F{Eit, 0, t rY(t,o,i,^)L^) = F{Tt (C, lo)) 

= Fi^,iu)^fiO,^,iu) 

and the property is proved when y — 0. 
For y e arbitrary, 

fiY{t,y,i,iv),E{t,y,^,iu),uj) = f{Y{t,0,^,Tyiu) + y,E{t,0,^,TyU;),iu) 

= fiY{t, 0, TyLU), E{t, 0, C, Tj^t^), TyLU) 

= /(0,C,r,^) = /(2/,^,c.) 

according to the result in the case y = 0. 



□ 



mvanance 



Remark 2.1. We mention here an important but easy consequence of the relations ^ and the 
of the measure P w.r.t. Ty, y g ; for any stationary function f = f{y,^,uj) — F{^,TyUj), F g 
L°°(M^ x^),we have 

E[f{Y{t, y, y, •), O] = E[F{T,{S„ ■))] 

for all t > 0, y, ^ E M.^ ; in particular, the left-hand side of the above equality does not depend on y. 
This property was used in the proof of lemma \2J\ 



Remark 2.2. Let us precise a little what happens when the function F e Lj^^(R^, L^(il)). In that case, 
F G L^(K^ X il; /ic) for all c > 0. Consequently, for any c > 0, we can define the function Fc associated 
to F by corollary \ 2. 1[ 

It is then easily proved that for any < c < c' , Fc = Fc' , ^c-almost everywhere. Setting An — 
{{^,Lo) e Supp^„; F„{y,() ^ F„+i{y,^)}, and A = U^^q^"; see that /ic(^) = for all c > 0. 
Moreover, for all (£,,uj) € x \ A, for all integers k,l such that {^,uj) G Supp/^^ n Supp^/, we have 
Fk{£,,uj) — Fi{£_,uj). We can thus define a function F{^,u}) on R^ x il \ A by 

F(^,a;) — Fn{£,,uj) for any n e N such that {^,uj) G Supp/i„ 

We then now that ^ 

F(Tt(^,c.))dt^F(C,c.) (9) 

as T ^ CO, and the convergence holds in L^{^c) for all c > 0, and fin almost everywhere for n G N. 
Eventually, setting 

1 r'^ 

B -.^ {{^,uj) eM.^ X n\A; — F{Tt{^, lu)) dt does not converge towards F{^, lu)} 

T Jo 

it is easily proved that f^dB) — for all c > (the equality is true for c e N, and is then deduced for 
c > arbitrary because the family of measures (fic) is increasing in c). 

Eventually, we have found a function F £ L\^JW^ ,L^{Vl)), independent of c, such that ([9]) holds i: 



m 



L^{lic) and fic-almost everywhere for all O 0. 



Remark 2.3. The construction above allows us to make more precise what we mean by projection P 



let f = f{y,^,u}) be a stationary function, f G L°° , Ll^{R^ , L^fl))) , and set F{^,uj) = /(0,C,w) £ 
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L]^^{R^ , L'^ (n)) . We can then associate to F a function F € L]^^{R^ , L'^ (n)) such that ^ holds in 
L^(^c) for all c (see remark W^) . We set 

P{f){y,^,uj) ■.^F{i,Tyu). 

It follows from corollary \2.1\ that P{f) is invariant by the hamiltonian flow ([U, and thus satisfies the 
constraint equation. From now on, we take this definition for the projection P, instead of the one given 
in the introduction. Notice that, for all y G and ^c- almost everywhere, 

P{f){y,£„uj) = lim i / F{Tt{tTyLo)) dt 

1 

= Jim Tf; / f {y{t,^-,^,TyUj),^{t,(^-,£.,TyUj),TyUj) dt 
= lim ^ f f{Y{t, y,^,u;),E{t,y,^, cu),cu) dt 
And we also give a more precise definition of S}{y,S^,ijj) : let 

Then ^ 

C»(2/,e,^) = ^(C)(y,e,^)= lim ^ I ~{t,y,^,u)dt 

a.e. and in L^{^c) for all < c < oo . 

Eventually, we mention here a property that will be used in the proof of the theorem; with the same 
notations as above, let 

<l^{T,y,£,,uj) = F {Tr{£„TyUj)) . 
Then is a solution of the evolution equation 

dr(t) + i- Vy<p - Vyu ■ Vj^ = 0, 

with initial data ^(r — 0, y, ^, lu) = f{y, £_,uj) = F Tyu). 

3 The general A^-dimensional case 

This section is devoted to the proof of theorem [TJ The proof is divided in three steps : first, we study 
the case of an initial data which does not depend on x, then the case when the initial data only depends 
on X (and not on y,^,uj), and eventually, we treat the general case. 

3.1 First case : /o does not depend on x 

Here, we assume that /o — fo{y,^,u}) G Ll^^{R^ ■,L°°{Ry x ft). Recall that /o is stationary, i.e. /o(2/ + 
z,£^,uj) — fo{y,£,,TzUj) a.s. in oj, for all {y,z,^) e R'^^. In the rest of the subsection, we set 

Fo{^,Lo) :=/o(0,C,c^) 

and ^ 

Fo{^,Lo) lim 1 / ^(C,^)) dt, fo{y,^,io) = Fo{^,TyU;). 

Notice that F e Ll^{R^ ; L°° (n)) , and thus F e L\R'^ x n; for all c> 0. 
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In that case, 



t X 



fo[Y [ -,-,Cuj ] -,-,^,uj] ,uj 



t X 



e e 
t 



e e 

/o(r(-,0,e,T^c^),S(-,0,C,r^t^),T^c^ 



T^U) - Fr 



rit,x,tuj) 



= /o(f'C'-) + {^o(T|(e, 

In accordance with theorem [J, we set 

<?(r, y, t u) = (Fo - Fo) (TV TyCu)) , 

and r'^ = 0. Then g satisfies the microscopic evolution equation ([3]) thanks to the remark at the end 
of the preceding section. Moreover, g{T) G K-^ by definition of K-*- and because P {Fq (T^- {^,TyUj))) = 

Fo{^,TyUj). 

It only remains to check that 

'"fix \ 
g -,— ,f,cj dt ^ ase^O 

in ii\,^(R^,Li(]R^ X for all T > and c> 0. 

The invariance of the measure P with respect to the group of transformations {Tx)xeW' (see remark 
12.11) entails that 



dfj.ci^,uj) 



1 



Foim,iu)) dt-Fo{^,iu) 



and the term above goes to as e ^ according to corollary 12.11 and is independent of a; G R^. There 
remains to check that /o = P{fo)- This follows directly from remark [231 Thus theorem [T] is proved in 
the case when /o does not depend on the macroscopic variable x. 

The following remark will prove to be useful when treating the general case : 

Remark 3.1. // /o G L°° , then for any function a G L°°{{0, oo) x x M.^ x Q), stationary in y, we 
have 



9 [ ~i ] dt as e 



in Lj^^iRx, L\R^ X n, n^)) for all T > and c> 0. 



Indeed, prove the property first for a = ai{t)a2(jj,^ji^), with ai, 02 G L°^\ For a arbitrary, take a 
sequence as with S > 0, converging to a in L\^^, and such that 



ns 

= ^af(i)a^(2;,e,t 

fc=0 



with af, a2 in L°° . The property is known for as, and it is thus easily deduced for a. 



3.2 Second case : /o = fo{x) 

Unlike the preceding subsection, we now focus on the case when /o only depends on the macroscopic 
variable x. In order to simplify the analysis, we assume that fo G W^'°°{R^) (the case when /o is not 
smooth in x will be treated in the next subsection). In that case, 

nt,x,t^) = fo(£Y(^-,^,t^ 
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Hence we have to investigate the behavior as e — > of 



We prove the following 

Lemma 3.1. Let T > arbitrary. As e vanishes, 

' t X 

e' e 



in L°°((0,T) X ;Li(Rf x n,fic)). 
Proof. Let us write, for t > 



eY{i,^,^,uA~x + t^^ (l'^' 



OJ 



e 
-t- 
t 



Let < a < r arbitrary. For a < i < T, we have 

/ eY(i^,i,u)-x + te(-,i,u) 

jR"xn f Jo 



< T sup 



dnc{£„i^) 

Li(R«xn,^c) 



and the upper-bound vanishes as e — > for any a > thanks to corollary 12.11 Notice that the upper- 
bound does not depend on x, hence the convergence holds in L°°{R^ ; L^{^c)). 
We now have to investigate what happens when t is close to 0; notice that 



sup 



lie" (g''?'^) llLi(K"xf2.Mc) - ^0 



where the constant Co only depends on N and c. Similarly, for alH > 0, 



sup 



Li(R"xn,/ie) 



< Co. 



Hence, if < i < a, we have 



sup 



dfj.c{^,LL!) < 2aCo. 



Eventually, 



eY 



L~((0,T)xR«;Li(/ie)) 



< inf <^ 2Coa + T sup 

0<Q<T T->a 



< 
1 

T 



and the lemma is proved. 



□ 
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We easily deduce that theorem [T] is true when /o € W^'°° {M.^) with 

f{t,x,y,£,,uj) := fo{x - if "(j/, f , w)), g = 0, 
r'(t, X, f , Lu) f{t, X, ^,Lu)- f (t, x,^,t 

and it is easily checked that / satisfies P(f) = f, f{t = 0) = P{fo) — fo (since /o is independent of y 
and f ) , and that / is a solution of the evolution equation (|4]) . 

3.3 Third case : fo arbitrary 

We now tackle the case of an arbitrary stationary function fo € Ll^^{M.^ x M^, L°°(]R^ x il)). We begin 
with the case when 

fo{x, y, uj) = a{x)b{y, f , uj), 

with a e W^i'°°(R^) and b e Lj^^ (Rf , (R^ x 17^^^ rc^mN 
directly from the two first subsections. Indeed, let 



i^^V^^^ , ^ xn))nL°^'{R^ xR^ xn),b stationary. This case follows 



and 



fit, x, y, t ^) = aix - tS}{y, i, lj)) P{b){y, t w), 



g{t, x- T, y, C, Lo) = a{x - tS}{y, ^ lo)) (b - P{b)) (Triy, f, lo)). 



It is already known that / and g satisfy ([4]), that f{t,x,-) G IK, and that g satisfies ^ thanks to the 
preceding subsections and the fact that ^'(j/,^,^) is invariant by the Hamiltonian fiow {Y,E). Notice 
that it is capital here that the coefficient ^"(2/, f , w) in the transport equation (0]) belongs to K. 

There remains to check that g{t, x; r, •) S IK-*-, that the remainder goes to strongly in Lj^^ and that 
g{t, x; t/e, x/e, uj) goes weakly to in the sense of theorem [TJ First, notice that a{x — t£}{y, f , lo)) S K 
and (6 - P{b)) {Tr{y,^,uj) e K^. Thus, a{x - t^^)P{b) = P{a{x - t^^)b) almost everywhere (because 
a{x — if'*(0, f , cj)) is invariant by the semi-group Tr), and consequently 

g{t, X- T, y, e, Lo) = [a{x ~ te)b - P {a{x - ti^)b)] {TA^, TyU;)). 

Hence g{t,x;T) E K-^ a.e. 
Then, setting 



t X 

e e 



r%t,x,£_,u;) = f{t,x,^,uj) - f (^t,x,^,£„u?j -g (^t, 
we have to prove that goes to strongly in Ll^^. We compute the difference 

f^t,x,£_,u;) - f (t,x, -^,S.,u?j - g (^t,x; ^, 



t X 



£ £ 



= a[eY[-,-,t^ b[Y [-,-,tu; ,^[-,-,C,LU ,u; 



t X 



£ £ 



t X 



£ £ 



-a(x-tS}{^-^,i,uj)) P{b)[-^,i,uj) 

-a (x - te« (f , u))[b- P{b)] ( J- f > , S ( J, ^, ^, c 



The right-hand side of the above equality is bounded by 



IHk^,~||&||i 



£Y 



and thus converges to as e — > in L°°((0, T) x R^ ; L (R^ x 51, /ic))) according to the second subsection. 
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Moreover, it is easily proved that as e 0, 







9 {t,x;-,-,^,uj] dt^O 



e e 



strongly in Ll^^{B.^ x M.^,L^{il)) thanks to remark ISTTl Hence theorem [T] is proved in that case. 

Now, let /o G Llci^^ X ,L°°(M^ x n)) arbitrary, and set Fq{x,^,lu) /o(x,0,^,w). Take a 
sequence of functions Fn £ L^(R^ x x Q) such that 

• -> i^o as n -> cx) in L^^iM.^ x x Q); 

• For all n e N, there exist functions e n W^'^{R^), hi £ n {R^ x n), 1 < k < n such 
that 

n 

F,,{x,^,cu)=J2ak{^)bm,^) a.e. 

Let be the solution of HI) with initial data _F„ (x, ^, r|w) , and let /„ = fn{t,x,y,S,,uj), gn = 
gn{t,x;T,y,^,uj) be the functions associated to by theorem [T] for all n. 
Let f{t, X, y, uj), g{t, x; r, y, ^, u) be the solutions of the system 



dt 



f 



Pif) = L Pig) = o, 



drg + £, ■ ^yg - V.yU(t/, uj) ■ V^.g = 0, 
fit = 0) = P{fo), git = 0, x; r = 0, ^, lu) = [/o - P(/o)] (x, y, c^). 

We have already proved that /„, gn satisfy the above system. We denote by Fq, Fn, the functions 
associated to Fq, F„ respectively by corollary 12.11 so that P(/o)(a;, y, ^, = Foix,^,TyUj), and /„(i = 

0,X,y,^,Oj) = Fnix,^,TyUj), gnit = 0,X,T = 0, ?/, ^, w) = (F„ - Fn)ix,£„TyUj). 

We use the following lemma, of which we postpone the proof : 

Lemma 3.2. Let g^ be a solution of ^ with initial data go {x, |, ^,a;), and go e L;^^(R^ xR^, L°°(R^ x 
ri)) stationary. Then for all R, R', T > 0, for all t e [0, T], 

\g%t,x,^,uj)\ dxd^< ||.go||Li(if^^^,,L~(R«xn)) 

where 



KtmM' = e X R^,|x| < i? + ry/i?'^ + 2u„,„ |^| < ^ R'^ + 2u, 

Similarly, if g is a solution of ^ with initial data go G ^[^^(R;^' x R^,L°°(R^ x fi)), i/ien 

\git,x,y,£,,Lv)\ dx < / ^ |5o(a;,2/,'?,w)| 

x<R J x<R+Ty/$^+2u^ 



Consequently, with Cu^t,^ ■— {x G R^, |a;| < R + + 2ujnax}5 we have 

\\fit,-,y,C^) - fnit,-,y,^,^)\\L^Ba) - W^oi' , t Ty^^) ~ ^ni' , ^, TyUj)\\ ^^^^^ 



< 



\Fo~Fn\i;tryiu) 



LHCR,T.i) 



||/(t,a;,y,C,w) - fnit, x,y,^,uj)\\^,^s^^T^N ^^^^^ ^^^^^^^^^ < \Fo - F„|(x,C,w) 
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In the last inequality, we have used property ([8]). 
And similarly, 

\\f{t,X,^,Uj) ~ fn{t,X,t^)\\L^(BnxBR,xn) < 11^0 " ^^n I iLHif^,,.^, ,L-(n)) . 

I \g(.t, X] T, y, ^, uj) - gn{t, x; r, y, ^, t^)|lLi(B„xRf xo,d:r d^,{i,c.)) < 2 | |Fo - F„(x, t^)! lLi(c„_^_^xRf xn,di.<iMe(«.c^)) • 



The above inequalities are true for all t e [0, T] and for all t > 0. 
Set 



t X 

£ £ 



Then for all t e [0, T], for all ti e N, setting c = iff 



1 E.'2 



< writ)- mm 

+ ~ 5"(0llL=°((0,oo)^xR^;)Li(Bi, xR|'xn,dx dMc(C,i^)) 

+\K{mLHB 

< 4 ||i^o - -^""1^1(0^ J, ^xRf xO,d2;dMc(C.'^)) + lkn(OllLi(BBXS^,xn) 

Thus r= ^ as £^ in L°°([0, oo); L]i„^(Mf x R^;L'^{n)). 

There only remains to check that g{t,x;t/e,x/e,^,uj) dt goes strongly to in L\^^ norm as £ 
vanishes; this result follows immediately from the same property for gn and the above inequalities. 
Therefore, we skip its proof. 

Proof of Lemma [37S[ First, let us recall that 



t X 



e £ 



t X 



£ £ 



t X 



£ £ 



and the .Jacobian of the change of variables 



t X 



£ £ 



£ £ 



is equal to 1. 

On the other hand, since 



we have 
and 

Thus 



<W\^\^ + 2u, 



\f%t,x,C,u;)\ dxdC 



xEBr, ^EBr, 



xeBn, ^EBr, 



t X 



e £ 



t X 



£ £ 



t X 



£ £ 



< 



< 



sup 

x€Br, JGB^, jjGR" 



t X 



£ £ 



t X 



£ £ 



dx d^ 



R"xR« 



l(sF(|,f,«,.),H(i,f,«,.))eif,,«,„, sup 



sup \fa{x,y,^,Lu)\ dx dS, 



t X 



£ £ 



dx 



„ f t X 



£ £ 



da; 
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The proof of the other inequahty goes along the same hnes. 



□ 



4 The integrable case 

In this section, we treat independently the periodic and the stationary ergodic case. Indeed, some results 
of the periodic case are no longer true in the stationary ergodic setting, and the results which do remain 
valid are not proved with the same tools. 

Let us make precise what we mean about "integrable case" : in the periodic case, we take a function 
u{y) which has the form 

N 

1=1 

where each function Ui is periodic with period 1 {1 < i < N). The Hamiltonian H{y,^) can be written 

1=1 

where Hi{yi,£) — jl^ip + Ui{yi) {1 < i < N). And the Hamiltonian system Q becomes 

i = (11) 

Yi{t^O)^y„ S,(i = 0)=6. 

Thus it is enough to investigate the behavior of each one-dimensional Hamiltonian system ifTTj) individ- 
ually, and for most calculations, we can assume without loss of generality that N = 1, and we drop all 
indices i. However, for the calculation of the projection P, a more thorough discussion will be needed, 
and we will come back to the case when TV > 1 in the corresponding paragraph. 

In the stationary ergodic setting, expression lfTO|) can be transposed in the following way : assume 
that fl = HfLi^li, where each ^li is a probability space, and assume that for 1 < z < A^, an ergodic group 
transformation, denoted by {Ti^y)y^M, acts on each fli. 

Then for u) = (wi, • • • ,u)n) e fl, and y = {yi, • • • , yw) e M^, we set TyUj := {n^y^uji, ■ ■ ■ ^TN^y^UN). 
And we assume that the function u can be written 

N 

u{y,i^) ^^Ui in^y^uji) , 

i=l 

where Ui £ L°°{ni) for all 1 < i < A^. The same remarks as in the periodic case can be made, and thus 
we will only consider the case iV = 1; note that in the stationary ergodic case, we are unable to compute 
the projection P when iV > 1. 



4.1 Periodic setting 

The goal of this subsection is to give another proof of the results of K. Hamdache and E. Frenod in [7], 
based on the study of the system 

y = -s, 

E = u'{Y), (12) 
Y{t^Q)^y, S - y e M, e e R. 

The Hamiltonian H{y,^) = ^ICP + u{y) is constant along the trajectories of the system (fT2l) . so that 

^\Eit,yO\^ +uiYit,y,0) ^ H{y,0- 
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We now fix y,^ G M^. Without any loss of generality, we assume y G [~\, 5), and we set £ :— 
H{y,£^). The above equation describes the movement of a single particle in a periodic potential u, with 
< M < Mmax- It is well-known that there are two kinds of behavior, depending on the value of the energy 
£ : if f < Mmax, the particle is "trapped" in a well of potential around y, and Y{t) remains bounded as 
t ^ 00. In that case, the trajectories in the phase space are closed curves. If £ > Wmax, the trajectory 
of the particle is unconstrained and \Y{t)\ ^ 00 as t ^ 00. We study more precisely these two cases 
and their consequences on the expression of the projection P in the following; we refer for instance to 
[2] for further calculations and results about Hamiltonian dynamics and ordinary differential equations 
in general. 

4.1.1 Expression of ^"(y, ^) 

We begin with the case when H(y,^) < Wmax- In that case, u{y) < < Wmax- By continuity of 

the potential u, there exists y- < y and y+ > y such that H{y,^) < u{y±) < Wmax, and the periodicity 
of u allows us to choose y± such that |y+ — y-| < 1- Then y_ < Y{t,y,£^) < y^ for all t > 0. Indeed, 
assume that there exists t > such that Y{t,y,^) > y+ > y — Y{t ~ 0,y,^). Since the trajectory Y is 
continuous in time, there exists < tg < t such that Y{t = tg, y, ^) = y_,_, which is absurd since 

H{Y{t = to, V. 0, S(t = to, V, 0) = H(y, > u{Y{t = to,y, 0) > H{y, 0- 
Thus Y{t,y,£^) is bounded. Since 



^«(y, = lim ^ I y, - - 1™ \; I Y{t, y, dt = lim : 

we deduce that £}{y,C) = for all y,^ such that H{y,^) < Umax- 
We now study the case H{y,^) > Umax- Since 

|r(t, y, = 2 {H{y, - <Y{t, y, C))) > 2(iJ(y, - Umax) > 
Y does not vanish for i > 0. Consequently, 



T 



S(t, y, i) = -Y{t, y, C) = sgn(e) ^2{H{y,i) - u{Y{t,y,i))l 



and since y,^) — y| > ^2 (iJ(y, ^) — Umax)^, \Y{t)\ ^ 00 as i ^ 00. We immediately deduce that 
S(t,y,<^) is periodic in time: indeed, there exists to > such that 

i"(io,y,0 =2/-sgn(0- 

And 



S(t = to, y, - sgn(e) v/2(i/(y,0-"(>^(io,2/,0)) = sgn(e) ^2 (if(y, - u(y)) = ^ = S(t = 0, y, , 
so that for s > 0, 



r(to + s,y,0 = l'(s,y,C)-sgn(^), 
S(to + s,y,^) = S(s,y^), 



and ^ is periodic with period to. 
Consequently, 



lim rE{t,y,Odt = ^ f"E{t,y,Odt. 

T^oc 1 Jq to Jo 



But 



["E{t,y,Odt = - f''Y{t,y,Odt 
Jo Jo 

= -(r(to,y,0-2/) 

= sgn(0. 
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Thus we only have to compute to- With this aim in view, we use the change of variables s = Y{t), with 
Jacobian ds — Ydt (recall that Y{t, y, = — sgn(^)-\/2 {H{y, ^) — u{Y{t, y, ^))) ), in the formula 



to 

dt 

y{to) 



Yit=o) -sgn(e)V2 iHiy,0-u{s)) 
/•a-sgn(C) Y 



ids 



v/2 iHiy,0-uis)) 
ds 



zds 



Eventually, we deduce 
where 



^2{H{y,0~uis)) 

e«(y,0 = sgn(0¥>(i?(y,0), 

1 



(£) = \/2l£>„„ 



We close this paragraph with a calculation which allows us to express in terms of the homogenized 
Hamiltonian H. The result we will obtain will be justified in more abstract and theoretical terms in the 
last subsection, using arguments similar to those of the theory of Aubry-Mather. 

First, let us recall the expression of the homogenized Hamiltonian H (see [10]) : we have 



H{y,^) = - \^\'^ + u{y), withinfM = 0, supu = Un 



and thus 

In other words, setting 

we have 
Hence, 

and 



ij if p < (^V2(u,„ax - u) 

2 I A if bl > ( v/2(umax -u)) , where \p\ = ( y^2{uma^ - u) + X 



[0,oo) 
X 



[0,oo) 



V20 



u) + X 



-ff(p) = "max + -l|p|>e(0)6' ^(bl)- 



H'{p) ^ sgn(p)-l|p|>e(o) — 



9'(A) 



1 



2 \ V2("max -U) + X/ ' 

0-i(b|) = 2(7?(p)-w) V|p| >0(O), 
IpI > 9(0) ^ H{p) > M„,ax Vp. 

Gathering all the terms, we are led to 

1 



H'{p) = sgn(p)\/2l5(p)>„_^ 
= sgn(p)(p {H{p)) 
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Thus, the final expression is 
where p is such that 

Hip) = H{y, ^) V Umax, sgn(p) = sgn(^). 



4.1.2 Expression of the projection P 

We also mention here how to find a general expression of the projection P in the special case iV = 1, 
and we explain how to generalize this expression in some particular cases when iV > 1. Recall that if 

/ = f{y,0 e Lloc(My X ^f) is periodic in y, then 







r 

T^oo T , 

and the limit holds almost everywhere and in L^([0, 1) x R^, nic), with dmdy, = '^H{y,^)<c dy d^. 

We begin with the case H{y,£^) > Mmax- We have seen in the previous paragraph that there exists 
to > 0, which depends only on H{y, ^) such that for all t > 0, for all fc G N 

Y{t + kto, y, = y{t, y, - fcsgn(e), E{t + kto, y, = S(t, y, 0- 
Thus f{Y{t),E{t)) is periodic in time with period to, and 

Pif) (y, = r /*" /(^(*' ^)) 

to Jo 

We use once again the change of variables s — Y{t), so that 
fiYit,y,0,^it,y,0)dt 



y 



f{s,sgniO^/2{H{y,0-uis))) ^ ds 

-sgn(0\/2 {H{y,£_)~u{s)) 



/(•,sgn(0V2(i?(y,0 -"(•)) 



1 



^2{H{y,0-u{-))/ 
Pif){y,0^fisgn{0,H{y,0) (13) 



And eventually, 
with 



V(£-"), 

We now focus on the case < £^ < Umax- In order to simplify the analysis we assume that E ^ 
{u(jj) ; u has a local extremum at y} (this set is finite or countable), and that 

Vy e K, u{y) = 0^u has a local extremum at y. 

In that case, it can be easily proved that Y{t, y, ^) is periodic in t; this follows directly from the fact 
that the trajectory in the phase space is closed (see [2J). Indeed, pushing a little further the analysis of 
the previous paragraph, we construct z± such that 

|z+ — z_| < 27r, z_ < z+, 
u{z±) = £, 
z- <y < z+, 
u{z) < £ Vze(z_,z+). 
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Then the particle starting from y with initial speed — ^ reaches either or z_ in finite time; the speed 
of the particle is at that moment since 

\Y\'^2{S~u{Y)), 

but its acceleration is —u'{z±) ^ 0, so the particle turns around and goes back in the reverse direction. 
It then reaches the other extremity of the interval (z-,z+) in finite time, and the same phenomena 
occurs. Hence after a finite time to, the particle is back at its starting point y with the same speed 
Consequently, the movement of the particle is periodic in time with period io- Thus, we have 

where ii > is arbitrary. It is convenient to choose for ti the first time when the particle hits Z-. In 
that case, it is easily seen that to is twice the time it takes to the particle to go from z_ to z+, so that 

^0 f*''''"^\, n 1 /-, 1 \ 



and 



Jt, \ J2(£-u) 



^2 [E-uY ^2 {£-u)^ 



ti+ta 



1 . 1 



f{Y{t, y, 0, S(t, y, 0) dt - ( l„<£/(s, )- ^ 



Gathering all the terms, we are led to 



lu<£ 



2(lu<£- ' 



(14) 



Expressions lfT3|) and (fT4|) are compatible with the ones in [7]. 

Let us now come back to the case when iV > 1, and take a function ip{y, ^) = ¥'1(2/1, Ci) ■■■ ¥'jv(2/jv, ■Cjv) 
where each (fi is periodic with period 1. We want to compute the limit 



1 '•^ 



TJo 



miYiit,yi,^i),:ii{t,yi,^i)) ■ • • .^at (1^ (t, yAr, ^at), ^jv(t, yjv, ^at)) dt. 



In general, knowing the behavior of each trajectory (1^, S^) independently is not enough to compute such 
a product. However, here, we recall that each function (pi{Yi{t, yi, £,i),'E,i{t, yi, ^i)) (1 < ? < N) is periodic 
in time. The period depends only on Hi{yi,^i) and on the function Ui. More precisely, setting 



/•I I 

Jo W£-Ui(z) 



\j£ - Ui{z) 

(pi{Yi{t,yi,^i),Ei{t,yi,^i)) is periodic in time with period Ti{Hi{yi,£,i)). 
We can thus use the following result : 

Lemma 4.1. Let /i, • • • , Jn & L°°{M.) such that fi is periodic with period di, and set (fi) — j- J^' fi. 
Assume that 

^ + ••• + ^^0 V(fci,...,fc^)eZ^\{0}. (15) 

fi tl]y 



Then as T ^ 00, 

^ ^ /i(t)---/^(i)rfi->(/i)---(/jv) 



TJo 
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Sketch of proof . By density, it is enough to prove the lemma for /i, •, /at S C°°(R). Write fi as a Fourier 
series (the series converges thanks to the regularity assumption), and use the fact that for all a ^ 0, 

T Jo 

□ 

In the present setting, we deduce the following result : 

Proposition 4.1. Let ip : (y,f) i-^ ",51(2/1,6) ' ' ■'PNiyN,^N), where G i^^(Rj^ x M«). 

Let (2/,C) e [0,1)^ X R^, and let 9i = e,{yi,i,) ^ T,{H,{yi,Q) for 1 < i < N. Assume that 
(6'i,--- ,9n) satisfy condition lfT5|) . Then 

P{ip)iy,0 ^ Piiipi)iyi,^i) ■ ■ ■ PNiipN)iyN,^N) (16) 
where each Pi is the projection in dimension 1 with potential Ui, given by expressions pTS)) and lfT4|) . 
In particular, when the set 

{(y,e) e [0,1]"^ X M^; (0i(yi,Ci),-- - ,^jv(2/jv,?jv)) satisfy condition dm)} 

has zero Lebesgue measure, equality (fT6|) holds almost everywhere. It can then be generalized to arbitrary 
functions Lp € L|^j.(M^ x R^) (always by linearity and density). The correct expression of the projection 
P is then 

P = PioP20---oPn, (17) 

where each projection Pi acts on the variables (jji,^i) only. Notice that all projections Pi thus commute 
with one another; hence the order in which they are taken is unimportant. 

We wish to emphasize that on the open set {(y, £_) € ]R^^,Vj S {!,••• , N} Hi{yi, 6) > max Wi}, the 
expression (fT7| is true. Indeed, the function Ti is strictly decreasing on (maxw^, +00), and thus the set 

{(£1, ■■■ ,£N)e R^; £^ > maxu, and 3fc e \ {0}, A:i/ri(£:i) + • • • + kN/TN{£N) = 0} 

is countable. As a consequence, the set 

{(y,^) e 'E?^ ,Hi[yi,S^i) > maxw^ Vi and (6'i (yi, Ci), • • • ,0NiyNAN)) satisfy condition (Ull)} 

has zero Lebesgue measure. 

However, let us mention here that in general, condition ifTSj) cannot be relaxed : indeed, assume for 
instance that Ui = uj :— u for i ^ j and assume that the function u is such that 

3yo > 0, u{y) = y^ for |y| < yo, 

and u{y) > y^ if y e [-i, i] \ [-yo,yo]- 
Then if |f | < y/yo, we have 

T{£) ^ / , dy^2 , dz =: Tq 

Thus, if Hi{yi,S^i) < y/yo, then {Yi,Ei){t,yi, £^i) is periodic with period Tq. Notice that Tq does not 
depend on the energy Hi{yi,^i) 

In that case, the function (f{Y{t), '^{t)) is also periodic with period Tq. Thus we have to compute the 
limit of ^ 

^j^ fl{t)---fN{t)dt 
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as r ^ oo, where the ft are arbitrary functions with period Tq. It is then easily proved that 

r 

T 



1 f 

Tf; I fi{i)- ■ ■ fN{t) dt ^ ^ ai,fci • • • aAr,fc„ (18) 

k e Z^, 

fci H h fcAT = 

where ^ 

-to Jo 

In general, the right-hand side of (fT8| differs from ai^ ■ ■ ■ aj^f^, and thus 

for (y, ^) in a neighbourhood of the origin. 
4.2 Stationary ergodic setting 

In the stationary ergodic setting, some of the expressions or properties above are no longer true. The 
most significant difference occurs when the energy H{y,£^) < Umax; indeed, in that case the particle is 
not necessarily trapped, depending on the profile of the potential u. Hence, in the rest of the subsection, 
we focus on the case Hiij,^) > Wmax- In that case, the movement of the particle is unbounded and has 
many similarities with the periodic case. In particular, the particle sees "all the potential" during its 
evolution, and this will be fundamental in the use of the ergodic theorem. 

4.2.1 Expression of ^"(y, ^, 

This paragraph is devoted to the proof of proposition [TTT] in the stationary ergodic setting. We refer for 
instance to [12] for conditions on the existence of correctors for all P G K in the case of a general coercive 
hamiltonian. In the present case, there exist correctors if 

{y G M; u{y, u) = Umax} 7^ 

a.s. in Lo G fl. 

Remark 4.1. We wish to point out that the expressions in the periodic and in the stationary ergodic case 
when H{y,^,uj) > Umax dre exactly the same (compare proposition \l.l\ and the end of paragraph \4-l.l^ . 
This expression, and more precisely, the equality = H'{P) for some P, is in fact strongly linked to 
Aubry-Mather theory. Indeed, 

1 r^u. Y{T,y,^,^)-y 



e«(y, ^, u) = lim ^ f Eit, y, lo) dt = - lim 



T^oo T 

and ^'^{y,^,uj) is thus (up to a multiplication by —1) the rotation number associated to the Hamiltonian 
flow starting at (y,C)- I'f^^ interested reader should compare our proposition \m\ to lemma 2.8 in [5] or 
theorem 4-1 in J6], and our proof to the ones in these articles. We refer to J^lEj for further references 
to Aubry-Mather theory and its applications to partial differential equations. 

Proof of proposition In all the proof, we fix y,£,,u! such that 77(y,^,w) > Umax, and we set P = 
P{y,^,u)). Let Q G M be arbitrary, and let u be a corrector, i.e. 

H{y,Q + Vyv[y,u;),Lo)^H{Q). 

Then according to the theory of viscosity solutions, for all (y, ^, w) G K x M x il, for all T > 0, 

v{y, iv) < viY{T, y, ^, u;)) + [ L{Y{t, y, ^, u;), E{t, y, ^, u;),u;)dt + Q [Y{T, y, ^, uj)-y]+ H{Q)T. (19) 
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Hence 



and 



L{Y{t,y,^,uo),S{t,y,^,oo),uj) dt> 



1 



viy)-viY{T,y,^,co)) YiT,y,^,u)-y 



-H{Q)T (20) 



-\Y{t, y, C, ^)| = H{y, ^,lu)-u> H{y, ^, lu) - u^^^ > Vi > 0. 
Thus there exist constants a,(3 > depending only on H{y,^,Lu) and Wmax, such that 

Y{T,y,^,u;)-y 



< a < 



T 



< /3 VT > 0. 



Consequently, Y{T) — > cx) as T — > oo and 

- v{YiT, y, ^, w)) _ ^(2;) - viY{T, y, u)) Y{T, y, - V 



T 

(remember ([6])). 

On the other hand, as T — > 00, 



r(T,2/,C,c.)-y 



as r ^ 00. 



Y{T,y,i,uj)-y _ 1 



T 



(21) 



(22) 



Hence, passing to the limit in l(20|) . we derive 

P(L)(y, ^, c^) > Qe"(2/, ^, ^) " H{Q) VQ G M. 

Thus 

P{L){y,^,u:)>L{^\y,i,oj)). 

In order to prove the proposition, we have to find a special Qo G IR such that equality holds in |(19 
This will entail that 

P(L)(y,e,c.) = sup {QS}{y,^,u:) - H{Q)) = L {£}{y,^,^)) , 

and the sup is obtained for S}{y,£^,uo) = H'{Qo). 

Let us thus prove that with Q = P = P{y,^,Lu), equality holds in ifTOl) . 
First, notice that 



is a viscosity solution of 
and as 2/ ^ 00 

By definition of H, 

consequently. 



v{y, w) sgn(P) ^2{H{P) - u{z, uj)) dz ~ Py 
H{y,P + VyV,Lo)=H{P), 



1 

y Jo 



2{H{P)-u{z,uj))dz^ E 



2{H{P) - U) 



2{H{P) - U) 



= \Pl 



E 



sgn(P) / d2{H{P) - u{z, uj)) dz ~ Py 



as y — > cx), a.s. in oj. Thus v satisfies ([6]), and v G L°°(ri;C^(R^)). Thus the method of characteristics, 
for instance, can be used to prove that equality holds in lfT9|) . with (y,5) replaced by any couple = 
(y',P + Vyv{y' ,ui)). We have to prove that we can take (y',C') = iv^O- First, notice that 



Wyv{y\Lu) = sgn(P)^2(i7(P) - u{y',uj)) - P, 

and thus sgn {P + \/yv{y' ,u} j) — sgn(^') = sgn(P) = sgn(^). Hence, take y' = y. Then — \£_' 
because H{y, ^) = H{P) = H{y,C) by definition of Thus ^ = and equality holds in 



l\2 



□ 
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4.3 Expression of the projection P 

The same method as in the periodic case can be used in order to find the expression of the projection P 
when H(y,^,io) =: £ > Wmax; indeed, in that case, remember that 

P(/)(y,e,a;)= lim ^ / J {Y{t,y,i,io),^{t,y,^,io),io) dt 

and we can use the change of variables 

dt=^dY= , ^ dY 

Y -sgn{OV^{£-u{Y,Lo)) 



in order to obtain 



f{Yit,y,tuj),E{t,y,^,uj),uj) dt 

= / /(^,sgn(0V2(^-u(2:,w)),a;) -—-j=====dz. 

Jy y ' -sgn(0V2(^-'«(z,w)) 

Since the group transformation {jx) is ergodic, and Y{T) ^ oo as T ^ oo, for all f > Wmax, 

— — — / f iz,sgn{^)^2{£ - u{z,ujj),iv] == dz ^ 

^{-l-J-yJv ^ ^ -sgn{Q^J2{£ -u{z,uj)) 

V y -sgn{i)^/2{£ - u{Si,u)) 

Thus, we obtain 

P{f) (y, ^, ^) = {y, ^)f{sgn{0, H{y, i,Lo)), 

where 

1 



{rjV^{£-uiO,uj)),u?j 



r]^/2{£-u{0,uj)) 
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